Mon. Not. R. Astron. Soc. 000, 1-?? (2002) Printed 2 February 2008 (MN BrgX style file vl.4) 



Modelling the infrared galaxy evolution using a 
phenomenological approach 



G. Lagache/ H. Dole,^ J.-L. Puget^ 

^ Institut d'Astrophysique Spatiale, Bat. 121, Universite Paris XI, 91405 Orsay Cedex, Prance 
'^Steward Observatory, University of Arizona, 933 N Cherry Ave, Tucson, AZ, 85721, USA 



Accepted 2002 August 14, Received 2002 June 6 



ABSTRACT 

To characterise the cosmological evolution of the sources contributing to the infrared 
extragalactic background, we have developped a phenomenological model that con- 
strains in a simple way the galaxy luminosity function evolution with the redshift, 
and fits all the existing source counts and redshift distributions. Cosmic Infrared Back- 
ground intensity and fluctuations observations, from the mid-infrared to the subniil- 
limetre range. The model is based on template spectra of starburst and normal galax- 
ies, and on the local infrared luminosity function. Although the Cosmic Infrared Back- 
ground can be modeled with very different luminosity functions as long as the radiation 
production with redshift is the right one, the number counts, and the anisotropics of 
the unresolved background, imply that the luminosity function must change dramat- 
ically with redshift, with a rapid evolution of the high-luminosity sources (L>3 10^^ 
Lq) from z—0 to z=l which then stay rather constant up to redshift 5. The derived 
evolution of the IR luminosity function may be linked to a bimodal star formation 
process, one associated with the quiescent and passive phase of the galaxy evolution 
and one associated with the starburst phase, trigerred by merging and interactions. 
The latter dominates the infrared and submillimetre ouput energy of the Universe. 
The model is intended as a convenient tool to plan further observations, as illustrated 
through predictions for Herschel, Planck and ALMA observations. Our model pre- 
dictions for given wavelengths, together with some usefuU routines, are available for 
general use. 

Key words: galaxies: evolution - infrared: galaxies - galaxies: starburst - galaxies: 
general - cosmology: observations 



1 INTRODUCTION 

The discovery of the Cosmic Infrared Background (CIB) 
(Puget et al. 1996; Fixsen at al. 1998; Hauser et al. 1998; 
Schlegel et al. 1998; Lagache et al. 1999; Lagache et al. 
2000; see Hauser & Dwek 2001 for a review), together with 
recent deep cosmological surveys in the infrared (IR) and 
submillimetre (submm) has opened new perspectives on our 
galaxy formation and evolution understanding. The surpris- 
ingly high amount of energy contained in the CIB showed 
that it is crucial to probe its contributing galaxies to un- 
derstand when and how the bulk of stars formed in the 
Universe. Thanks to ISO (Kessler et al., 1996) - mainly 
at 15 ^lm with ISOCAM (Cesarsky et al. 1996), and 90 
and 170 /xm with ISOPHOT (Lemke et al. 1996) - and 
ground-based instruments - SCUBA (Holland et al. 1998) 
and MAMBO (Bertoldi et al. 2000) at 850 and 1300 ^m 
respectively - deep cosmological surveys have been carried 
out. It is thus now possible, to various degrees, to resolve the 



CIB into discrete sources (e.g. Kawara et al. 1998; Barger 
et al. 1999; Elbaz et al. 1999; CariUi et al. 2000; Juvela et 
al. 2000; Linden- Vornle et al. 2000; Matsuhara et al. 2000; 
Bertoldi et al. 2001; Dole et al. 2001; Elbaz et al. 2002; Scott 
et al. 2002). The striking result of these surveys concerns 
the evolution of the IR and submm galaxy population. The 
source counts are high when compared to no, or moderate, 
evolution models* for IR galaxies (Guiderdoni et al. 1998; 
Franceschini et al. 1998). Therefore, it has been necessary 
to develop new models in the IR. Very recently, several em- 
pirical approaches have been proposed to model the high 
rate of evolution of IR galaxies (e.g Devriendt & Guider- 
doni 2000; Wang & Biermann 2000; Gharry & Elbaz 2001; 
Franceschini et al. 2001; Malkan & Stecker 2001; Pearson 
2001; Rowan-Robinson 2001; Takeuchi et al. 2001; Xu et al. 



* 'No-evolution': the co-moving luminosity function remains 
equal to the local one at all redshift 
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2001; Balland et al. 2002; Wang 2002) which fit all existing 
source counts, redshift distribution and CIB intensity and 
fluctuations, although often not all of them. We present in 
this paper a new model whose preliminary results were pub- 
lished by Dole ct al. (2000). The originality of this model 
was to empirically separate the evolution of the starburst 
galaxies with respect to the normal galaxies, the current ob- 
servations implying a strong evolution of the bright part of 
the Luminosity Function (LF). It was shown for the first 
time that only the starburst part should evolve very rapidly 
between z=0 and z=2, the evolution rate being much higher 
in the IR than in any other wavelength domain. 

We present here a more sophisticated and detailed ver- 
sion of the first model (Dole et al. 2000). Our philosophy 
is to build the simplest model of the LF evolution, easily 
deliverable, with the lowest number of parametres but ac- 
counting for all observationnal data. We stress out the point 
that we include the CIB fluctuations levels as measured by 
ISOPHOT (Lagache & Puget 2000; Matsuhara et al. 2000) 
and IRAS (Miville-Deschenes et al. 2002) as an extra con- 
straint. Recent observations strongly suggest that the bulk 
of the optical and IR cxtragalactic background is made of 
two distinct galaxy populations (see Sect. 2.2). Therefore, we 
restrict our model in the wavelength domain 10 - 1500 nm, 
our goal being to quantify the evolution of IR galaxies. The 
paper is organised as follows: we first summarise our present 
knowledge on the evolution of IR galaxies, and on the nature 
of the sources contributing to the Extragalactic Background 
(Sect. 2). Then, we present the ingredients of the model 
(Sect. 3). In Sect. 4, we discuss the galaxy templates used in 
the model. We then present the parametrisation of the local 
LF (Sect. 5). In Sect. 6 arc given the results of the model 
(evolution of the LF, evolution of the luminosity density, 
number counts, z-distribution, CIB intensity and fluctua- 
tions). And finally the model is used for predictions for the 
Herschel And Planck surveys (Sect. 7.2 and 7.3 respectively), 
gives requirements for future large deep survey experiments 
(Sect. 7.4) and predictions for ALMA observations (Sect. 
7.5). A summary is given in Sect. 8. 



2 OUR PRESENT KNOWLEDGE 

2.1 The strong evolution of IR galcLxies: 
observationnal evidences 

There have been, in the last few years, strong observationnal 
evidences indicating extremely high rates of evolution for IR 

galaxies. 

First, galaxy evolution can be observed through its im- 
print on the far-IR EB (Extragalactic Background). Weakly 
constrained even as recently as 6 years ago, various observa- 
tions now measure or give upper/lower limits on the back- 
ground from the UV to the millimetre (mm) waveband (e.g. 
Dwek et al. 1998; Gispert et al. 2000; Hauser & Dwek 2001). 
The data show the existence of a minimum between 3 and 
10 fiin separating direct stellar radiation from the IR part 
due to radiation re-emitted by dust. This re-emitted dust 
radiation contains at least a comparable integrated power 
as the optical/near-IR, and perhaps as much as 2.5 times 



more. This ratio is much larger than what is measured lo- 
cally (~30%). The CIB is thus likely to be dominated by a 
population of strongly evolving redshifted IR galaxies. Since 
the long wavelength spectrum of the background is signifi- 
cantly flatter than the spectrum of local star-forming galax- 
ies, it strongly constrains the far-IR radiation production 
rate history (Gispert et al. 2000). The energy density must 
increase by a factor larger than 10 between the present time 
and a redshift ~ 1-2 and then stay rather constant at higher 
redshift (till ~ 3), contrary to the ultraviolet radiation pro- 
duction rate which decreases rapidly. 

Second, several deep cosmological surveys at 15, 90, 170, 
850 and 1300 /im have resolved a fraction of the CIB into dis- 
crete sources. For all surveys, number counts indicate a very 
strong cosmological evolution of IR galaxies, in total power 
radiated but also in the shape of the LF. This is particularly 
obvious at submm wavelengths where the EB is dominated 
by high luminosity galaxies (see the SCUBA and MAMBO 
results) . The high rates of evolution exceed those measured 
in other wavelength domain as well as those observed for 
quasars and Active Galactic Nucleis (AGNs) . 

Finally, high rates of evolution are suggested by the de- 
tection of the CIB Poissonian fluctuations at a high level at 
60 and 100 /um with IRAS (Miville-Dechenes et al., 2002) 
and 170 /xm with ISOPHOT (Lagache & Puget 2000; Mat- 
suhara et al. 2000). For example, Matsuhara et al. (2000) 
give the constraints on the galaxy number counts down to 
35 mjy at 90 fj,m and 60 mjy at 170 fmi, which indicate the 
existence of a strong evolution down to these fluxes in the 
counts. 



2.2 Sources making the Extragalactic 
Background 

2.2.1 Optical versus IR and submm EB sources 

Recent observations show that the bulk of the optical and 
IR EB is made by two distinct galaxy populations (e.g. 

Ausscl ct al. 1999). Therefore, one of the key question is 
wether the dusty star forming galaxies are recognizable from 
optical/near-IR data alone. 

In the local Universe, Sanders & Mirabel (1996) show 
that the bolomctric luminosity of IR galaxies is uncorre- 
lated with optical spectra. The color excess derived from the 
Balmer lines ratio does not significantly depend on the IR 
luminosity, IR color or optical spectral type (Veilleux et al. 
1999). In fact, neither the moderate strength of the heavily 
extinguished starburst emission lines, nor their optical col- 
ors can distinguish them from galaxies with more modest 
rates of star formation (Elbaz et al. 1999; Trentham et al. 
1999; Poggianti et al. 1999). However, recently, Poggianti & 
Wu (2000) and Poggianti et al. (2001) show that the inci- 
dence of e(a)^ sources in the different IR selected samples 
seems to suggest that the e(a) signature might be capable of 
identifying from optical data alone a population of heavily 

^ Galaxies with strong Balmer absorption lines and [Oil] in emis- 
sion (Poggianti & Wu (2000) and references therein) 
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extinct starburst galaxies. Reproducing the e(a) spectrum 
requires the youngest stellar generations to be significantly 
more extinguished by dust than older stellar populations, 
and implies a strong ongoing star formation activity at a 
level higher than in quiescent spirals. 

At intermediate redshift (0.3 < z < 1), one of the 
main information on IR galaxies comes from identifications 
of ISOCAMdeep fields. For example, Flores et al. (1999a,b) 

presented results of a deep survey of one of the CFRS fields 
at 6.75 and 15 /xm. At 15 /xm, most (71%) of the sources 
with optical spectroscopy are classified as e(a) galaxies. The 
FIR luminosities of the Florcs ct al. c(a) galaxies arc be- 
tween 5.7 10"'^" and 2 lO'^'^ L©. This is the first confirmation 
of the IR luminous nature of e(a) galaxies. More recently, 
Rigopoulou et al. (2000) using VLT spectroscopy found that 
optical ISOGAM counterparts in the HDF-south are indis- 
tinguishable from the dusty luminous e(a) galaxies. 

All these studies at low and intermediate redshift seem 
to show that IR galaxies predominantly exhibit e(a) signar 
tures in their optical spectra. On the contrary, in optically 
selected surveys of field galaxies, o(a) spectra are present 
but seems to be scarse (less than 10%, see for example Pog- 
gianti et al. 1999). In fact evidence for high e(a) incidences 
are found in merging or interacting systems or active com- 
pact groups. A complete study of the IR emission of e(a) 
galaxies is still to be done. 

At much higher redshift, there are other evidences that 
the optically selected samples and bright IR samples (e.g. 
the SCUBA blank field sample) are different. For example. 
Chapman ct al. (2000) have performed submm photometry 
for a sample of Lyman Break galaxies whose UV properties 
imply high star formation rates. They found that the inte- 
grated signal from their Lyman break sample is undetected 
in the submm. This implies that the population of Lyman 
break galaxies does not constitute a large part of the de- 
tected blank field bright submm sources. 

In conclusion, it is clear that the optical and IR EBs are 
not dominated by the same population of galaxies. There- 
fore, we restrict our model in the wavelength domain 10- 
2000 /^m, our goal being to quantify the evolution of the IR 
galaxies. 

2.2.2 The GIB sources 

In the HDF-N field (HST Deep Field North), Ausscl et al. 
(1999) and Elbaz et al. (1999) find that 30 to 50 % of ISO- 
GAM galaxies are associated with optical sources showing 
complexe structures and morphological pecularities. More- 
over, Cohen et al. (2000) show, in the HDF-N, that more 
than 90% of the faint ISOGAM soxaccs are members of con- 
centrations. This shows that past or present interactions or 
merging play a large role in trigerring the IR emission of 
galaxies. All the studies of ISOGAM field sources show that 
the bulk of the CIB at 15 ixm comes from galax;ies which 
have bolometric luminosities of about 10^^-10^^ L©, high 
masses (~ 10^^ Mq) and redshift between 0.5 and 2. They 
experience an intense stellar formation (100 Mo/yr) which 
appears to be uncorrelated with the faint blue galaxy pop- 



ulation dominating the optical counts at z~0.7 (Ellis 1997; 
Elbaz et al. 1999). In galaxy clusters all /SO CAM sources 
are found preferably at the periphery where there is still 
some star formation (Biviano et al. 1999). 

At longer wavelengths, source identifications are much 
more diflScult. Thus, characterising the nature of the galaxies 
is a long time process. However, we already have some in- 
dications. The FIRBAGK survey at 170 /im detected about 
200 galaxies (Dole ct al. 2001) making less than 10% of the 
CIB. Schematically FIRBAGK sources comprise 2 popula- 
tions: one cold and nearby (L~ 10^ - lO'^^ Lq) and one 
cold or warm very luminous (L~ 10^^ L©) with redshift 
lower than 1.2. The optical spectroscopy of the brightest 
FIRBAGK sources reveals an 'IRAS-likc' starburst nature 
(Dennefeld et al. in prep; Patris et al. 2002) with a moder- 
ate star-formation rate (10 Mq/jt). These results are very 
similar to those of Kakazu et al. (2002) who found that 62% 
of the Lockman Hole 170 jj,m sources are at rodshit below 
0.3 with luminosities lower than lO^'^ Lq (based on Arp220 
spectrum), the rest being mostly ultraluminous IR galaxies 
with redshift between 0.3 and 1. In their sample, the 170 
^m sources appear also to be powered primarily by star for- 
mation. 

In the submm, the main indication comes from the SGUBA 
deep surveys (e.g. Hughes et al. 1998; Barger et al. 1999; 
Eales et al. 1999; Scott et al. 2002). These surveys suggest 
that faint 850 fim sources are mostly ultraluminous galax- 
ies at typical redshift between 1 and 4 (e.g. Eales et al. 
2000). SGUBA sources above 3 mjy account for 20-30 % of 
the EB at 850 yum. The present data show that the bulk of 
the submm EB is likely to reside in sources with 850 pm 
fiuxcs near 0.5 mJy. Barger ct al. (1999) estimate that the 
FIR luminosity of a characteristic 1 mJy source is in the 
range 4-5 10^^ Lq, which is also the typical luminosities of 
sources making the bulk of the CIB at 15 ^m. Moreover, as 
for the 15 /xm sources, several groups have suggested that 
the submm sources are associated with merger events (e.g. 
Small et al. 1998; Lilly et al. 1999). All these results show 
that the integrated power of the LF at redshift greater than 
~0.5 must be dominated by sources with luminosities of a 
few 10^^ Lq while the local LF is dominated by sources with 
luminosities of the order of 5 10^° Lq . 

2.2.3 The AGN contribution to the IR ouput energy 

At low-z, Veilleux et al. (1999) and Lutz et al.(1999) show 
that most of the IR sources are powered by starbursts. The 

AGN contribution appears dominant only at very high lumi- 
nosities (L>2 10^^ Lq). Also using FIRBA GK nearby bright 
source spectroscopy, we find that less than 10% of sources 
show AGN signs (Dennefeld et al., in prep; Patris et al. 
2002). The same conclusions are reached by Kakazu et al. 
(2002). At intermediate z, from optical and X-rays studies 
of ISOGAM sources making the bulk of the CIB, several 
groups show that 15 /im sources are mostly starburst galax- 
ies (e.g. Fadda et al. 2002). At much higher z, the main in- 
dication comes from X-ray observations of SGUBA sources 
(Barger ct al. 2001; Fabian ct al. 2000; Hornschemeier et 
al. 2000; Severgnini et al. 2000). All these observations of 
submm galaxies in the X-rays are consistent with starburst- 
dominated emission. However, recently. Page et al. (2002) 
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present a result of SCUBA observations of eight X-ray ab- 
sorbed active galactic nuclei from z=l to z=2.8 and find, for 
half of them, a 850 jj,m. submm counterpart. Nevertheless, 
the high 850 /xm fluxes (greater than 5.9 mjy) suggest that 
these sources are liyporlurninous galaxies. Such galaxies do 
not dominate the IR output of the Universe. 

Considering the whole CIB energy budget, and based on 
the assumptions that 10% of the mass accreting into black 
hole is turned into energy and that the black hole masses 
meEisured in the HDF (Ford et al. 1998) are typical of galax- 
ies, the AGN background energy would be at the order of 
10% of that from stars (Ealcs ct al. 1999). These calcula- 
tions arc highly uncertain but arc supported by the work of 
Almaini et al. (1999). And recently, more direct evidences 
have been obtained. For example, Severgnini et al. (2000) 
show that the 2-10 keV sources making at least 75% of the 
X-ray background in this band contribute for less than 7% 
to the submm background. 

On the modelisation side, several groups (e.g. Xu et 
al. 2001; Rowan-Robinson 2001) show that the starburst 

galaxy evolves much more rapidly than the AGN-dominatcd 
sources, making the AGN contribution to the CIB relatively 
small. 

In conclusion, it is clear that AGNs do not dominate the 
IR output energy of the Universe. Therefore, we will con- 
sider hereafter that IR galaxies are mostly powered by star 
formation and we use Spectral Energy Distributions (SED) 
typical of these star- forming galaxies. The differences in SED 
for the small fraction of AGN dominated IR galaxies would 
not affect significantly the results of our model which is built 
with only two galaxy populations ('normal' and starburst) 
defined by their SEDs. 



The differential counts at a given flux S and wavelength 
A=Ao write as: 

dN f f dN . dz^^ . ,^ ^ 

We then obtain the integral counts: 

N{> S) = j ^dS (5) 

the CIB intensity produced by all sources with S<Smax'- 

Ic.s = I S-dS (6) 

and the CIB intensity fluctuations (the shot noise) from 
sources below a given detection limit So (that could corre- 
spond to either the confusion or the sensitivity limit) mear 
sured by the level of the white noise power spectrum: 

Pf^^^ = j\^^dS JyVsr (7) 

Wc build the simplest model with the lowest rmmbcr of 
parametres and ingredients which fits all the observations. 
Wc first fix the cosmology (f2A=0.7, Oo=0.3 and h=0.65) 
from the combination of the most recent CMB anisotropy 
measurements (e.g. de Bernardis et al. 2002) the distance- 
luminosity relation of type la supernovae (e.g. Perlmutter et 
al. 1999) and the measure of galaxy distances using Cepheids 
(e.g. Preedman et al. 2001). We then construct the 'normal' 
and starburst template spectra : at each population and lu- 
minosity is associated one spectrum (see Sect. 4). We fi- 
nally search for the best evolution of the LF (Eq. 3) that 
reproduces the number counts (Eq. 4-5), the CIB and its 
fluctuations (Eq. 6 and 7 respectively), assuming that the 
LF evolution is represented by the independent evolution of 
the two populations. A rather remarquable result, as will be 
seen later, is that two populations only can fit all the data. 



3 INGREDIENTS OF THE MODEL 

At a given wavelength A=Ao, the fiux 5*1^ of a source at 
redshift z, as a function of the rest-frame luminosity L,^ (in 
W/Hz) can be written: 



S'.(L,2,A = Ao) = 



il + z)x K{L,z) X L,{L,\ = \o) 



(1) 



where Di is the luminosity distance and K(L,z) is the K- 
correction factor defined as : 



K{L,z) = 



i^(z=0) 



(2) 



This correction is specific of the spectrum of the considered 
population at given L and z. In practice, the rest-frame lu- 
minosity Li, is convolved by the band pass filter centred on 
A=Ao. The number of sources per solid angle and redshift 
interval is: 



^^^(L,.) = A.o(L,.)x(l + .)3x- 



(3) 



where dV/dz is the differential volume clement fixed by the 
cosmology and No is the number of sources per unit volume 
and luminosity interval as a function of redshift. No is given 
by the LF. 



4 BUILDING GALAXY TEMPLATES 
4.1 Starburst galaxies 

The luminous IR starburst galaxies emit more than 95% of 
their energy in the far-IR. Spectra of such galaxies have been 
modeled by Maffei (1994), using the Desert et al. (1990) dust 
emission model and the observationnal correlation for the 
IRAS Bright Galaxy Sample of the flux ratio 12/60, 25/60 
and 60/100 with the IR luminosity (Soifer & Neugeubauer, 
1991). We start from this model and modify it slightly to 
better take into account recent observationnal constraints. 
The significant improvements are the following: 

• First we replace the Desert et al. (1990) PAHs template 
by the Dale et al. (2001) one, keeping the same amount of 
energy in the mid-IR. 

• We slighly modify the spectral shape in the near and 
mid-IR: we increase the PAH and Very Small Grains pro- 
portions (by a factor ~2) but add some extinction, slowly 
increasing with the luminosity. 

• Finally, wc broaden the far-IR peak (also a continuous 
change with the luminosity) and flattens slightly the long 
wavelength spectrum (Dunne al. 2000; Dunne & Eales 2001; 
Klaas et al. 2001). 
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Figure 1. Color diagrams from the model (continuous lines), 
compared with Soifcr & Ncugeubauer 1991 (diamonds) and 
Dunne et al. 2000 (stars) data points. 

Fig. 1 shows the (12/25, 60/100) and (60/850, 60/100) color 
diagrams compared with Soifer & Neugeubauer (1991) and 
Dunne et al. (2000) observations. We see a very good agree- 
ment between the templates and observations. The average 
luminosity spectra sequence is shown on Fig. 2. For the 'nor- 
mal' starbursting galaxies (Li_r <10^^ L©), we also check 
that the 7/15 versus 60/100 diagram was in good agree- 
ment with Dale et al. (2001). 

Such a representation of galaxy SEDs that assigns only 
one spectrum per luminosity does not take into account the 
dispersion of the colors (e.g. the 60/100 ratio) observed for 
a given luminosity (as in Xu et al. 2001 or Chapman et al. 
2002a). However, we do not have enough measured colors 
to do a statistical analysis of their variations with the lumi- 
nosities. The only colors measured for a large sample are the 
IRAS colors (principally the 60/100). Using only the 60/100 
color as a tracer of the dispersion for a given luminosity, 
may not improve the representation as for example two very 
different long-wavelength spectra can have the same 60/100 
color (e.g. NGC 7821 and NGC 0549, Stickel et al. 2000). 

4.2 Normal spiral galaxies 

For the 'normal' galaxies (i.e. standard IR counterparts of 
spiral galaxies with more than half of their energy ouput 
in the optical) , we take an unique spectrum derived mainly 
from the /^OPTOT serendipity survey (Stickel et al. 2000) 
and the nearby FIRBACK galaxy SEDs, together with 
longer wavelength data from Dunne et al. (2000) and Dunne 
& Eales (2001). 

The ISOPHOT serendipity survey heis revealed a pop- 
ulation of nearby cold galaxies (Stickel et al. 1998, 2000), 
under-represented in the 60 jjbva IRAS sample. The distri- 
bution of Ii7o /Iioo flux ratio shows that about half of the 
galaxies have a flux ratio between 1 and 1.5, indicating that 




10 1 00 1 000 
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Figure 2. Starburst model spectra for different luminosities: 
L=3. 10l2 Lq (continuous line), L=5 lO^l Le (dotted line), 
L=10ii L0 (dashed line) and L=10i" Lq (dotted-dashed line) 



the FIR spectra are mostly flat between 100 and 200 /xm. 
Very few show a downward trend in this wavelength range, 
this trend being typical of warm starburst galaxies (see Fig. 
2). Most important is the large fraction of sources (more 
than 40 %) which have I170/I100 >1.5, indicating a rising 
spectrum beyond 100 /um similar to that seen for example 
in the Milky Way galactic ridge, a property known for more 
than 20 years from early balloon-borne measurements (Serra 
et al. 1978, Silverberg et al. 1978). 

In the FIRBACK Nl field (Dennefeld et al., in prep), the 
brightest sources show mean IR color similar to that of 
Stickel et al. (2000): I170/I100 ~ 1.3 and leo/Iioo ~0.47. 
These objects are often associated with bright optical spiral 
galaxies (e.g. Fig. 3) with also typical 15 (im. fluxes such as 
I170/I15 ~42. Finally recent observations at 450 /xm (Dunne 
& Eales 2000) reveal also the presence, in normal galaxies, 
of a colder component than previously thought. 
Therefore, we take for the normal galaxie the "cold" tem- 
plate presented in Fig 4. This template have mean IR color of 
Ii7o/Iioo= 1.42, l6o/Iioo= 0.35, Ii7o/Ii5=45 and Iioo/l850= 
77. For the Mid-IR part of the template, we use the spec- 
tral signature that applies to the majority of star-forming 
galaxy presented in Helou (2000) . More observations around 
the maximum intensity (~ 100-200 (im) and in the submm 
are needed if we want to refine the template and describe 
its variations with the luminosity. 

Note that we do not make any evolution of the 'normal' 
and starburst template spectra with the redshift. As shown 
in Chapman et al. (2002a), the available data for high red- 
shift, far-lR galaxies do not reveal evidence for any strong 
evolution in the characteristic temperature of the color dis- 
tribution over 0< z <3. We can also test this hypothesis 
of 'no-evolution' using the two high-redshift sources N 1-040 
and Nl-064 detected in the FIRBACK Nl field (Chapman 
et al. 2002b). For these two sources, we do a blind search for 
both the luminosity and redshift using our template spectra, 
based on a test. For Nl-064 the 'normal' and starburst 
template gives the same with Log(L)=12.2 and z=0.5 and 
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Figure 3. Example of association of one bright FIRBACK 170 
//m source with a "cold" spiral galaxy (the optical image is 
from the DSS; the white circle corresponds to the 170 //m er- 
ror position). IR color ratios for this source are l6o/Iioo=0.6, 
Ii7o/Iioo=l-4 and Ii7o/Iib=51 
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Figure 4. Template spectrum for the normal galaxies (continu- 
ous line), compared with the template spectrum for the starburst 
galaxies (dash-dot line), for the same luminosity L=5 10^" Lq 



Log(L) = 12.8 and z=1.05 respectively. The starburst tem- 
plate gives results in good agreement with Chapman et al. 
(2002b). For Nl-040, the best is obtained for the 'normal' 
template with Log(L) = 12.2 and z=0.45, which is in per- 
fect agreement with Chapman ct al. (2002b). The template 
spectra give photometric redshifts in good agreement with 
the spectroscopic redshifts, suggesting no strong evolution 
of the galaxy SEDs over Q < z <1. 
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Figure 5. Luminosity function at 60 /im (model: normal galaxies: 
dot-dash line, starburst galaxies: dash line, total: continuous line) 
compared with Saunders et al. 1990 (squares). 
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Figure 6. Bolometric luminosity function at z=0 (normal galax- 
ies: dot-dash, starburst galaxies: dash and total: continuous) 



5 PARAMETRISATION OF THE LOCAL 
LUMINOSITY FUNCTION 

A detailed comparison of the luminosity function of IR 

bright galaxies with other classes of extragalactic objects 
has been done in Sanders & Mirabel (1996). The most strik- 
ing results are (i) the high luminosity tail of the IR galaxy 
LF is clearly in excess of what expected from a Shechter 
function and (ii) at luminosities below ^lO^'^ Lq, the ma- 
jority of optically selected objects arc relatively weak far-IR 
emitters (as shown by IRAS, ISOPHOT and SCUBA ob- 
servations). Accordingly, we decompose the IR LF in two 
parts: (1) the 'low-luminosity' part (dominated by normal 
galaxies) that follows the shape of the optical LF and (2) 
the 'high-luminosity' part dominated by starburst galaxies 
(note however that each part covers the whole range of lu- 
minosities). 
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We start from the local LF at 60 /im from Saunders et 
al. (1990) (Fig. 5). We convert the 60 /xm LF in a bolomet- 
ric LF (L=l-1000 /xm) using our template spectra: at the 
'low-luminosity' and 'high-luminosity' part, we assign the 
'normal' and starburst template spectra respectively^ . The 
bolometric LF (1-1000 /xm) used in the model is shown in 
Fig. 6. 

The parametrisation and redshift evolution of the two 
parts of the bolometric LF is done using the following: 

• For the normal galaxies: 

- A LF at z=0 parametrised by an exponential with 
a cutoff in luminosity Lento ff- 

^normaliL) = ^bol{L) X exp ( — ) (8) 

\ i^cntoff J 

- A weak number evolution 

• For the starburst galaxies: 

- A LF described by: 

^SB{L,z) = ^boi{L = 2.h lO^Lo) (9) 

^ 2-5 10"Le j 1—) 

- An evolution of the luminosity of the knee (Lknee) 
with the redshift 

- A constant slope of the LF at high luminosities 
(SBsiope), not redshift dependent 

- An evolution of the luminosity density with the red- 
shift, ip{z), that drives the LF evolution: 

$sb(L, z) = $sb(L, z) X (^g^) (10) 

- An integrated energy at z=0 as observed: 

SBnorm='p{z = 0) 



6 RESULTS AND DISCUSSIONS 

6.1 LF evolution derived from observations 

Although the number of parametres is quite low, it is too 
much time consuming to do a blind search through the whole 

paramctre space. We therefore search the best solution for 
the parametres near values expected from direct observa- 
tionnal evidences: 

• We fix the normal galaxy evolution (passive evolution) 
such as it nearly follows the number density evolution of 
optical counts up to z=0.4: 

^norrnal{L, z) = ^normal{L, Z = 0) X (1 + z) . 

We arbitrarily stop the evolution at z=0.4 and keep this 
population constant up to z=5 and then let the population 
decreases up to z=8 (see Fig. 8) 

t We have also checked the validity of the long-wavelength part 
of our template spectra by comparing the 850 /im LF with Dunne 
et al. (2000) 



• An estimate of Lknee is given by deep surveys at 15 and 
850 ^m: the bulk of the CIB at 15 /im and 850 /im is made 
by galaxies with L~l-5 10^^ L© (e.g. Barger et al. 1999; 
Elbaz et al. 2002) 

• We have indications on the evolution of the luminosity 
density, ip{z), for the starburst galaxies, directly inferred 
from the CIB spectrum shape by Gispert et al. (2000) which 
is used as a starting point to adjust the model. 

• We take SB 

norm — 10^ Lo/Mpc'^ which is roughly the 
value observed locally. 

• We take SBsiope= -2.2 (Kim & Sanders, 1998) 

The best evolution of the LF that reproduces IR number 
counts, redshift distributions and CIB observations is shown 
on Fig. 7. It is obtained with: 

• Lcutoff = 5 10^^ Lq 

• Lfc,iee(z) = 8 10^° X (1-l-z)^ L0 up to z=1.5 and then 

Lfcnee(z)= Lfcnee (z= 1 .5) . 

Moreover, to avoid a 'break' in the evolved luminosity func- 
tion (as in Dole et al. 2000), we modify the low- luminosity 
part of the starburst LF with the redshift (up to z=5) ac- 
cording to: 

^sb{L,z) = $sb(L = ) X ^ (11) 

where Lmax is the luminosity corresponding to the maxi- 
mum of $sb(L, «). 

We see a very high rate of evolution of the starburst 
part which peaks at z~0.7 (Fig. 8) and then remains nearly 
constant up to z=4 (as shown for example by Charry & 
Elbaz, 2001). Wo compare on Fig. 8 the co-moving luminos- 
ity density distribution as derived from the model with the 
Gispert et al. (2000) one. There is a good overall agreement. 
However, the model is systematically lower than the Gispert 
et al. (2000) determination for redshifts between 0.5 and 2. 
This comes from the fact that the CIB values at 100 and 
140 /xm used in Gispert et al. (2000) were slightly overes- 
timated (Renault et al. 2001) leading in an overestimate of 
the luminosity density distribution at low redshift. 



6.2 Comparison model-observations 

6.2.1 The number counts 

Fig. 9 shows the comparison of the number counts at 15 /xm 
(Elbaz et al. 1999), 60 /um (Hacking & Houck 1987; Gre- 

gorich et al. 1995; Bortin et al. 1997; Lonsdale et al. 1990; 
Saunders et al. 1990; Rowan- Robinson ct al. 1990), 170 /im 
(Dole et al. 2001) and 850 ^m (Small ct al. 1997; Hughes et 
al. 1998; Barger et al. 1999; Blain et al. 1999; Borys et al. 
2002; Scott et al. 2002; Webb et al. 2002) with the observa- 
tions. We have a very good overall agroomont (we also agree 
with the 90 /im number counts of Serjeant et al. 2001 and 
Linden- Vornle et al. 2001). 
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6.2.2 The redshift distributions: Need for a normal 'cold' 
population 

On Fig. 10 is shown the redshift distribution of resolved 
sources at 15, 60, 170 and 850 /im. The 15 /xm redshift dis- 
tribution is in very good agreement with that observed by 
Floros ot al. (1999b) and Aussol ot al. (1999). At 850 fim, we 
predict that most of the detected sources are at z > 2.5. At 
170 ^m, wc predict that about 62% of sources with fluxes 
S >180 mjy (4a in Dole ct al. 2001) are at redshift below 
0.25, the rest being mostly at redshift between 0.8 and 1.2. 
We know from the F/_R_SACi<' observations that the redshift 
distribution predicted by the model is very close to that ob- 
served at 170 fj,m: it is clear from these observations that 
we have a bi-modal z-distribution (Sajina et al. 2002; Den- 
nefeld et al. in prep). Moreover, very recently, Kakazu et 
al. (2002) published first results from optical spectroscopy 
of 170 /xm Lockman Hole sources. They find that 62% of 
sources are at z<0.3 with IR luminosities (derived using 
Arp220 SED) lower than 10^^ L©, the rest being at redshift 
between 0.3 and 1, which is in very good agreement with 
the model prediction. It is very important to note that the 
agreement between the model and observations can only be 
obtained with the local 'cold' population. A model contain- 
ing only starburst-like template spectra and 'warm' normal 
galaxies overpredicts by a large factor the peak at z~l (as 
in Gharry & Elbaz 2001 for example). 

6.2.3 The GIB and its anisotropics 

The predicted GIB intensity at specific wavelengths, to- 
gether with the comparison with present observations are 
presented in Table 2 and in Fig. 11. We have a very good 
agreement with the estimates at 60 (Miville-Deschenes 
et al. 2002), 100 jum (Renault et al. 2001), and 170 urn (Kiss 
ct al. 2001) and the FIRAS determinations (Fixsen et al. 
1998; Lagache et al. 2000). We are also in good agreement 
with the lower limit derived from 15 /itm counts (Elbaz et al. 
2002), combined with the upper limit deduced from high en- 
ergy gamma-ray emission of the active galactic nucleus Mkn 
501 (Renault et al. 2001). According to the model, sources 
above 1 mJy at 850 /im contribute for about 30% of the GIB. 
At 15 /xm, with the deepest ISOCAM observations (Altieri 
et al. 1999; Aussel et al. 1999; Metcalfe 2000) about 70 % of 
the GIB has been resolved into individual sources. 

Finally, we compare the model predictions for the GIB 
fluctuations with the present observations at 60, 90, 100 and 
170 fxm (Table 3). For the comparison we need to remove 
the contribution of the brightest sources that make the bulk 
of the fluctuations (sources with fluxes S > Smax). When 
Smax is quite high (at the order of 500 mJy - 2 Jy) , the fluc- 
tuations are dominated by the strongest sources, making 
the result very dependent on the accuracy of the evaluation 
of Smax- This is why the comparison observations/model is 
very difficult^. For Smax around 50-150 mJy, the fluctuar 
tions are dominated by the faint and numerous sources that 
dominate the GIB and the values do not depend critically on 

§ We cannot compajre the model predictions with the Kiss et al. 
(2001) results since we have no information on Smax 



the exact value of Smax- We see from Table 3 that although 
for some observations the model can be lower or greater by 
factor 1.5 in amplitude, we have an overall very good agree- 
ment. We stress out that reproducing the GIB fluctuations 
gives strong constraints on the LF evolution. For example, 
evolution such that we better reproduce the 850 /xm counts 
gives too high GIB fluctuations. Future observations with 
better accuracy will show if these minor discrepancies dis- 
appear or are indicative that the bright submm counts are 
overestimated due to a high fraction of gravitationally lensed 
sources (Perrotta et al. 2002), or are simply indicative that 
the present phenomenological model is too simple ! The level 
of the predicted GIB fluctuations for dedicated experiments 
(as for example Boomerang or Maxima), with respect to the 
cirrus confusion noise and instrumental noise will be dis- 
cussed in detail in Fiat et al. (in prep). 

At 170 ^m, it is clear from Fig. 12 that the redshift 
distributions of sources that are making the GIB and those 
making the bulk of the fluctuations are similar. The fluctu- 
ations axe not dominated by bright sources just below the 
detection threshold but by numerous sources at higher red- 
shift. Thus, in this case, studying the GIB fluctuations gives 
strong constraints on the GIB source population. This is true 
in the whole submm-mm range. The population of galaxies, 
where the GIB peaks, will not be accessible by direct detec- 
tion in the coming years. For example, SIRTF will resolve 
about 20% of the background at 160 ^m (Dole et al. 2002), 
PACS about 50 % at 170 urn (Sect. 7.2.2) and SPIRE less 
than 10% at 350 /xm (Sect. 7.2.1). 

In conclusion, we have seen that our model gives num- 
ber counts, redshift distributions, GIB intensity and fluctu- 
ations that reproduce all the present observations. It can 
be now used for future experiment predictions, in particular 
for Herschel, Planck and ALMA observations. For SIRTF, 
a complete and more detailed study, including sirrmlations 
and a detailed discussion on the confusion, is done in Dole 
et al. (2002). 



7 PREDICTIONS FOR FUTURE 
EXPERIMENTS 

7.1 The confusion noise from extragalactic sources 

The confusion noise^ is usually defined as the fiuctuations 
of the background sky brightness below which sources can- 
not be detected individually. These fluctuations are caused 
by intrinsically discrete extragalactic sources. In the far-IR, 
submm and mm wavelengths, due to the limited size of the 
telescopes compared to the wavelength, the confusion noise 
play an important role in the total noise bugdet. In fact, 
the confusion noise is often greater than the instrumental 
noise, and is thus limiting severely the surveys depth. The 
total variance of a measurement within a beam due to 
extragalactic sources with fluxes less than Sum is given by: 

^ We only consider the confusion noise due to extragalactic 
sources since, in the high galactic latitude cosmological fields, 
the cirrus confusion noise is negligible. 
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a' = y" f(e,4>)ded4>j^ '^^ s^^ds (12) 

where f{9, (/>) is the two-dimensional beam profile (in stera- 
dians), S the flux in Jy, and ^ the diflerential number 
counts in Jy"'^ sr"'^. We call Sum the confusion limit. 

The confusion noise can be determined using two cri- 
teria: the so-called photometric and source density criteria 

(see Dole et al. 2002 for a full description). The photometric 
criterion is related to the quality of the photometry of de- 
tected sources, the flux measured near Sum being severely 
affected by the presence of fainter sources in the beam. It is 
defined by the implicit equation: 

Slim = Qphot X Cr{Slim) (13) 

where qphot measures the photometric accuracy and is usu- 
ally taken between 3 and 5. 

The source density criterion is related to the complet- 
ness of detected sources above Sum directly related to the 
probability to lose sources too close to each other to be 
separated. There is a threshold above which the density of 
sources above Sum is such that a significant fraction of the 
sources is lost (it is impossible to separate the individual 
sources anymore). For a given source density (with a Pois- 
sonian distribution) N(>S), corresponding to a fiux limit 
Slim, the probability P to have the nearest source with fiux 
greater than Sum located closer than the distance dmin is; 

P{< em^n) = 1 - e-"'^^-" (14) 

dmin is the distance below which sources cannot be sepa- 
rated and is a function of the beam profile. If we note 9fw, 
the FWHM of the beam profile, Omin can conveniently be ex- 
pressed in unity of ^fvk, Omin = kx 6fw- As an illustration, 
simulatious of source extraction at the MIPS/SIRTF wave- 
lengths show that k=0.8 should be achievable (Dole et al. 
2002). Therefore, in the following, we fix k=0.8. We choose 
the maximum acceptable probability of not being able to 
separate the nearest source "Pmax = 0.1. In this case, the 
source density is equal to 1 / 17.30 (Table 1) and the cor- 
responding most probable distance is about I.IxOfw- 
Using Eq. 14, we can derive N(>S) and thus find the cor- 
responding Slim- Then, using Eq. 12, we compute cr. The 
source density criterion leads to a value of an equivalent 
<ldenaity = If qdensity is greater than standard values 

of qphot (3 to 5), then the confusion noise is given by the 
source density criterion. If not, then the photometric crite- 
rion has to be used to derive the confusion noise. The clas- 
sical confusion limit of 1 source per 30 beams corresponds 
to k=l and P^O.l. Nevertheless it can still lead to mediocre 
photometry for very steep logN-logS. 

The transition between the photometric and source den- 
sity criterion is around 200 /im, depending on telescope di- 
ameters. For example for SIRTF, at 24 and 70 /im, the source 
density criterion gives qdensity ~7 (Dole et al. 2002). In this 
case, the source density criterion has to be used to derive 
the source confusion and it leads to a very good photomet- 
ric quality. On the contrary, at the Planck wavelengths, the 
source density criterion gives qdensity ~1, leading the source 
detection limited by the photometric quality. Furthermore, 



Table 1. Number of beams per sources for different probability 
P to have the nearest source located closer than the distance 
6min = kx dpw with O'^^y ~ n/1.1 (valid for both Gaussian 
and Airy disk beams). 





P = 0.05 


P = 0.1 


P = 0.15 


k=0.8 


35.6 


17.3 


11.2 


k=0.9 


45.1 


22.0 


14.2 


k=1.0 


55.7 


27.1 


17.6 



this illustrates the limits of doing high signal-to-noise obser- 
vations 'to beat the confusion'. 

In the following, we derive the confusion noise for fu- 
ture long wavelength dedicated surveys. We assume that the 
sources are randomly distributed on the sky. The effect of 
the clustering on the confusion noise will be investigated by 
Blaizot et al. (in prep) using the hybrid model described in 
Guiderdoni et al. (2001) and Hatton et al. (2002). The num- 
ber counts at long wavelengths derived from the model and 
usefuU for the following predictions are shown on Fig. 13. 

7.2 The Herschel deep surveys 

The 'Herschel Space Observatory' (Pilbratt, 2001) is the 
fourth cornerstone mission in the European Space Agency 
science programme. It will perform imaging photometry and 
spectroscopy in the far-IR and submm part of the spectrum, 
covering approximately the 60-670 /im range. Herschel will 
carry a 3.5 metre diameter passively cooled telescope. Three 
instruments share the focal plane: two cameras/medium res- 
olution spectrometres, PACS and SPIRE, and a very high 
resolution heterodyne spectrometre, HIFI. 

In Table 4 are shown the confusion noises and limits at 
PACS and SPIRE wavelengths using the photometric and 
source density criteria. At the PA C5 wavelengths, the source 
density criterion leads to confusion limits 2 to 6 times higher 
than the photometric criterion with qphot=5. For SPIRE, at 
250 ytim, the limits obtained using the two criteria are very 
close and then, at longer wavelengths, the photometric qual- 
ity constrains the confusion limit. 

Extragalactic surveys will be conducted by both the 
PACS and SPIRE instruments. Up to now, the trade-off be- 
tween large-area, shallow versus small-area, deep/ultradeep 
observations has not been finalised. In the following two sec- 
tions are examples of extragalactic surveys that could be 
done. 

7. 2. 1 The SPIRE surveys 

For large-area scan-mapping observations, the current esti- 
mates of time needed to map 1 square degree to an instru- 
mental noise level of 3 m,]y (lainst) is about 1.7, 2 and 2.1 
days at 250, 350 and 550 /xm respectively (Matt Griffin, pri- 
vate communication). Three kinds of surveys could be done 
with SPIRE (Table 5): 
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(i) A very large-area (2 approximatly 14° x 14°, i.e. ~ 400 
square degrees) survey at the noise level of 5cri„st=100 mjy 
at 350 fj,m. (92 mJy at 250 /um and 102 mJy at 550 /um). 

(ii) A confusion-limited survey of about 100 square de- 
grees. To reach a Bainat noise limit of 14.1, 22.4 and 17.8 
mJy, one needs 192, 90 and 149 days at 250, 350 and 550 jjm 
respectively. 

(iii) A very deep survey, down to the confusion limit to 
extract as much information as possible about the underly- 
ing population. For example, mapping 8 square degrees to 
5ai„st= 7.5 mJy at 350 /im would take 64 days. 

We concentrate, in the following, on the 350 /im band (al- 
though SPIRE is observing simultaneously the three bands). 
The number of detected sources for each survey is given in 
Table 5 and the redshift distribution in Fig. 14. 

• The very large-area survey 

This survey is well suited to the Planck sensitivity (Table 
7). It would provide better positions for the Planck point 
sources and the combination with the Planck data would 
improve the spectral and spatial characterisation of fore- 
grounds. It will take about 18 days, which is a small amount 
of time comparing to the enormous progress it would bring 
for the component separation problem for all observations 
in this wavelength range. The very-large-area survey will de- 
tect mostly nearby sources but hundreds of object could be 
detected at z>l. It will also explore the cirrus component. 
Moreover, thanks to its surface coverage, this survey will 
detect the rare and high-z objects for which little is known 
today. 

• The confusion-limited survey: 

For the confusion-limited survey, there are mainly two peaks 
in the redshift distribution: one at z~0 due to the cold 
sources, the other between 0.8 and 3 being due to the star- 
burst galaxies. The immber of starburst galaxies detected 
between 0.8 and .3 doesn't vary much. Fig. 15 predicts the 
number of starburst galaxies that the confusion-limited sur- 
vey will detect at 350 jim per log-interval of luminosities 
(this figure does not include the cold population which con- 
tributes mostly locally). The evolution of the 10^^-10^^ Lq 
galaxies will be measured from 0.5 to 2.5. Galaxies with 
L~3 lO^'^ Lq will be only accessible at z~0.1. The nature 
and number of the highest luminosity objects are important 
to test cosmological theories, and thus populating the high 
redshift bins with enough objects to ensure <14% Poisson 
noise (corresponding to 50 sources) is a critical driver for 
the size of the sample, and accordingly, the area required 
for the survey. Based on our model, we see from Fig. 15 that 
to meet this ^^14% goal for galaxies in the range 10^^ - 10^'^ 
Lq at z<2.5, an area of 100 square degrees is the minimum 
required. 

• The very-deep survey 

The confusion-limited surveys will detect only ~7% of the 
GIB (1% for the vey-large area survey). In this case, study- 
ing the fluctuations will bring informations on the underly- 
ing source population. The redshift distribution of sources 
making the bulk of the GIB and the fluctuations at 350 yum 
are very similar. Thus, the fluctuations will be the unique 
opportunity to get information, particularly on sources that 



dominate the GIB (sources with L~3 10^^ L©), at redshifts 
where they cannot be detected individually. To detect and 
study them, a field with high signal-to-noise ratio is needed. 
Moreover, the size of the field should be large enough to try 
to detect the source clustering and not only be limited to 
the Poissonian noise detection (which is about 4300 Jy^/sr). 
Since the source clustering is expected at the scales between 
1 and 5 degrees (Knox et al. 2001), a field of about 8 square 
degrees is the minimum required. This field could be part of 
the confusion-limited survey which will be certainly split in 
several smaller area surveys. 

7.2.2 The PACS surveys 

For the whole field of view of ~1.75' x 3.5', the current 
estimates of time needed to reach 5c7inst=3 mJy is 1 hour 
(Albrecht Poglitsch, private communication). At 75, 110 and 
170 /im, the confusion limits are about 0.13, 0.89, and 7.08 
mJy, with cidensUy of 8.9, 8.7 and 7.1 respectively. To reach 
the confusion limit for one field of view, i.e 5(Ti„3t=0.1, 0.9 
and 7.1, we need 567, 11 and 0.18 hours at 75, 110 and 170 
fxm respectively. Since for PACS the confusion limits and 
the time to reach those sensitivities arc very different at the 
three wavelengths, three kinds of surveys could be done that, 
schematically, will probe the GIB in the three bands: 

(i) A shallow survey of 20 square degrees down to the 
confusion limit at 170 yum. This survey is dedicated to probe 
the GIB at 170 /um and study the correlations. Such a survey 
will take 88 days. 

(ii) A deep survey down to the confusion limit at 110 
lira. A good compromise between the covered surface and 
the time needed is a field of 25'x25' that will take about 67 
days. 

(iii) An ultra-deep survey down to the confusion limit at 
75 yum (or a little bit below the confusion limit). To map a 
5'x5' field, 96 days are needed. 

These three surveys correspond to about the same amount 
of time as the SPIRE surveys. Obviously, the Pj4CyS' surveys 
have to be done within the same areas as the SPIRE ones. 

PACS is observing simulaneously the 170 and 110 /jm 
channels or the 170 and 75 /xm channels. Ideally a combinar- 
tion 75/110 /um and 75/170 um would have been preferable 
since the 170 /um observation in the deep and ultra-deep sur- 
veys will not bring new science compared to the large area 
survey. On the contrary, observing the shortest wavelengths 
in the shallow survey will detect the very luminous, hot and 
rare galaxies that may be missed in the deep and ultra-deep 
surveys. 

The three surveys will detect thousand of sources at 
z~l (Fig. 16) and will probe most of the GIB source pop- 
ulation (they will resolve about 49, 77 and 87% of the GIB 
at 75, 110 and 170 ^m respectively, see Table 6). 
At 170 pm, the shallow survey will give an unpreccdent mea- 
surement of the evolution of the lO'^'^-lO'^'^ L© galaxies from 
z~ 0.25 to 1 and the evolution of the 10^^-10^^ Lq galax- 
ies from z~ 0.5 to 3 (with enough objects to ensure <10% 
Poisson noise) . Since half of the background is resolved into 
discrete sources at 170 um, one complementary approach is 
to reduce the surface of the shallow survey to have a better 
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signal-to-noisc ratio and thus study the underlying popu- 
lation. However, to study the correlation in the IR back- 
ground, a minimum of 8 Sq. Deg. is required. A survey at 
Smin = 8 mjy, corresponding to 5i7i„at=S.7 mjy, would take 
around 128 days for 8 Sq. Deg. Such a survey would give loss 
statistics for the resolved sources but would help in under- 
standing the whole CIB population. With a 8 Square Deg. 
field, the measure of the high luminosity source evolution 
would still be possible with a high degree of accuracy. 

In confusion, PACS will resolve about 80% of the CIB 
around 100 jjxa {SIRTF wiW resolve at most around 55% of 
the CIB at 70 fim and 20% at 160 ^lm, Dole et al. 2002). It 
will definitely resolve the question of the population mak- 
ing the CIB near its emission's peak. It will measure with 
unprecedented accuracy the history of the IR-traced star 
formation up to z~1.5. For the higher redshifts, the informa- 
tions will come mainly from the SPIRE surveys. Although 
SPIRE will resolve less than 10% of the CIB in the submm, 
it will bring an unprecedented information on the evolution 
of galajcies up to z~3 and also on the underlying population 
that will be very hard to detect from the ground due to the 
small area of the present surveys. 

7.3 The Planck all-sky survey 

For the Planck wavelengths, the confusion limits are given 

by the photometric criterionll . With qphot=5, they are about 

447, 200, 79.4, 22.4 and 11.2 mJy at 350, 550, 850, 1380 and 
2097 /im respectively (for detection only, qphot=S is better 
and loads to the following confusion limits: 251, 112, 44.7, 
14.1 and 6.31 at 350, 550, 850, 1380 and 2097 fim respec- 
tively). The confusion limit (with qphot=5) is above the 5cr 
instrumental noise at 350 /ini, comparable at 550 /im and 
then below at longer wavelengths. It can be compared to 
the values given in the 'HFI proposal for the Planck mission 
(1998)', Table 2.1. In this Table, the confusion limits have 
been computed using Q. rather than J (6, 4>)d9d<p in Eq. 
12. This leads to a systematic overestimates of Oconf by a 
factor 1.33. Correcting for this factor, the ratio of those es- 
timates to the present ones increases from 350 to 2097 /^m 
from 1.2 to 3. 

To compute the rmmber of sources expected in the 
Planck survey, we use a cut in flux equal to hatot such as: 

where aconf for Planck is given by Geo 

is an additonal noise due to unresolved sources with fluxes 

between Sam and hcinsf- 

f{e,^)ded4, J s^^ds (16) 

This additional term is only present when the 5ainst is 
greater than the confusion limit (as for example for the 
SPIRE veTy-\aige survey). 



II The density criterion leads to qdensity from 1.5 to 0.8 from 350 
to 2097 iJ,m respectively. 



Final sensitivities on point sources, together with the 
number of detected sources for Planck are given in Table 
7. At all wavelengths, the sensitivity of the survey is in the 
euclidian part of the number counts. Planck will not be able 
to constrain, with the resolved sources, the evolution of the 
submm galEixies but it will give an absolute calibration of 
the bright number counts, that will not be provided by any 
other planned instruments. Moreover, by covering the whole 
sky, it will probably detect the most spectacular dusty ob- 
ject of the observable Universe, as the hyperluminous or the 
strongly lensed starburst or AGN galaxies, £is well as ex- 
treme sources not included in the model. 
Note however that the sensitivities have been computed with 
an instrumental noise derived from a mean integration time. 
With the Planck scanning strategy, some high latitudes re- 
gions (where the cirrus contamination is low) will be sur- 
veyed more deeply, leading to an instrumental noise about 
3 times lower. In such high redundancy parts of the sky, 
Planck survey will be limited by the confusion noise (except 
at 2097 jim where the instrumental and confusion noises 
are on the same order). In those regions, Planck will pro- 
duce unique maps of the CIB fluctuations. CIB anisotropies 
arc mainly contributed by moderate to high redshift star- 
forming galaxies, whose clustering properties and evolution- 
ary histories are currently unknown. Planck observations 
will thus complement the future far-IR and submm tele- 
scopes from ground and space that will perform deep sur- 
veys over small area. These surveys will resolve a substantial 
fraction of the CIB but will probably not investigate the 
clustering of the submm galaxies since it requires surveys 
over much larger areas. 



7.4 Requirements for future experiments for very 
large deep surveys in the IR/submm/mm 
domain 

Since the far-IR and submm sources (1) are often associated 
with mergers or interacting galaxies and (2) their energy 
ouput is dominated by high luminosity sources at high red- 
shift (which might be related to the optical where the out- 
put is dominated by high brightness sources at high redshift 
(Lanzetta et al. 2002)), they have to be studied in detail in 
the long wavelength range and at high redshift (z>3) as a 
tool to understand the merging process and the physics of 
the flrst non linear structures. 

Two kinds of requirements have to be meet by the post 

SIRTF, Herschel and Planck surveys: 

(i) They have to find the interesting high redshift sources 
which are early mergers made of building blocks not affected 

much yet by star formation and evolution. The number of 
such sources have to be large enough to do statistical studies. 

(ii) Once these sources are found, the future experiments 
have to have enough sensitivity and angular resolution to 
study them in details. 

To quantify the first requirement, we compute the sur- 
face needed to detect more than 100 sources with luminosi- 
ties of 3 lO'^'^ and 3 lO'^^ L© in each redshift range (Table 
8). At high redshift, we need surveys of about hundred of 
square degrees to find enough high luminous objects to do 
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statistical studies. If enough area is covered, we need more- 
over a high sensitivity, for example we have to reach 0.21 
mjy at 850 fxm for L=3 10^^ L© galaxies (dominating the 
LF at z~l) and 1.9 mJy at 850 /im (or 0.9 mJy at 1300 ^m) 
for L=3 10^^ Lq galaxies (dominating the LF at z>2). The 
fluxes of typical 3 10^^ Lq galaxies at high-z are just, for the 
single-antenna telescopes, at the confusion limit. With the 
future wide-field imaging instruments on these telescopes, 
for example SCUBA-2 (Holland et al. 2001) and BOLO- 
CAM (Glenn et al. 1998), hundred of square degrees at the 
1 mJy level at 1300 /xm could be mapped in a reasonable 
amount of time. However, with 10/15 arsec beam, it will be 
very difficult to make optical identifications and follow-up 
observations at other wavelengths. 

To reach the sensitivity of the 3 10^^ L© galaxies at high-z, 
we need to have a very good angular resolution not to be 
limited by the confusion. Table 9 gives the angular resolu- 
tion, together with the telescope diameter needed to reach 
10, 30, 60 and 80 % of the GIB at 350, 850 and 1300 ^m. To 
resolve 80% of the background at 1300 nm, we need a tele- 
scope diameter of about 173 m (!) and a diameter of about 
113 m (!) at 850 /xm and 23 m at 350 /xm. In conclusion, 
finding the objects that are making the bulk of the GIB 
at long wavelengths will be a very challenging task! This 
leaves to an open question: how we find these objects? In 
the mid-IR, the Next Generation Space Telescope (NGST) 
will be a great tool. However, NGST observations will have 
two limitations: (1) the redshift up to which the dust com- 
ponent can be observed is limited to 4-5 and (2) the stellar 
component can be observed at higher z but the experience 
of combined optical/near-IR and far-IR observations shows 
how it is difficult to identify the galaxies which have most of 
their output energy in the Far-IR from optical and ncar-IR 
data alone. Therefore, the alternative today is to make the 
interferometres efficient enough to carry out large surveys. 

For the second requirement, i.e study in detail the 
physics of the objects, observations have to have enough 
sensitivity to observe sub-components of the merging ob- 
jects at high z (a 3 10^° sub-component at z=5 has 
a flux of about 0.014 mJy at 350 /im and of about 0.028 
mJy at 850 pm). Moreover, an angular resolution of about 
0.2 to 1 arcsec, together with spectroscopic capabilities are 
needed. For that, the submm and mm interferometres are 
the only tools, as shown by the recent spectacular obser- 
vations of higli-rodshift quasar with the IRAM intorfer- 
ometre (Gox et al. 2002). Lower luminosity sources will 
require the ALMA (http://www.eso.org/projects/alma/, 
http://www.mma.nrao.edu/) or SPEGS/SPIRIT (Leisawitz 
et al. 2001) interferometres. 

7.5 The case of ALMA 

Thanks to the negative K-corroction, high redshift sources 
are accessible from the ground at submm and mm wave- 
lengths. ALMA, a synthesis radio telescope (about 64 12- 
metre diameter telescopes) that will operate at submm and 
mm wavelengths will image the Universe with unprece- 
dented sensitivity and angular resolution from the high- 
altitude Llano de Ghajnantor, in northern Ghile. It will be 
one of the largest ground-based astronomy project of the 
next decade after VLT/VLTI, and, together with the NGST, 



one of the two major new facilities for world astronomy com- 
ing into operation by the end of the next decade. ALMA, 
with its angular resolution, great sensitivity and spectro- 
scopic capabilities will reveal in detail, in the high-z galaxies, 
the astrophysical processes at work. Moreover, ALMA will 
be free of limitation due to source confusion and will there- 
fore allow very faint galaxies to be detected. In this section, 
we discuss mostly ALMA abilities to find large enough sam- 
ples of high redshift sources to do statistical studies and 
probe the GIB source population. Of course ALMA will be 
also used to study their structure and physics. 

At 1300 ^m, to find enough ~ 3 10^^ Lq high-z sources, 
we need to cover a surface of at least 5 square degrees at a 
5cr level of about 0.1 mJy (Table 8). Such a survey will re- 
solve about 50% of the GIB. To resolve -80% of the GIB, 
one needs to reach a 5(t level of about 0.02 mJy (Table 9). 
Therefore two kinds of surveys could be considered: a large- 
area (~ 5 Sq. deg.) and an ultra deep (~ 10 arcmin^) survey. 
For both surveys, the compact configuration has enough res- 
olution not to be limited by the confusion. At 1300 /xm, a 5(t 
detection of 2.3 mJy is reached in 1 sec for a beam area of 
0.16 arcmin^ {'ALMA proposal for phase 2' and Blain 2001). 
This gives for the two types of surveys: 

• The large-area survey: 5 Sq. Deg., 5(Ji3oo = 0.1 mJy 
(50% of the GIB) 

A 1 square degree field requires 22 500 pointings, each with 
529 seconds of observations, for a total of 138 days. For 5 

Sq. Deg., 690 days are needed, i.e. 1.9 years. 

• The ultra-deep survey: 10 arcmin^, 5ai3oo = 0.02 mJy 
(80% of the GIB) 

A 10 arcmin^ field requires 625 pointings, each with 13225 
seconds of integration, for a total of 96 days. 

In conclusion, if we want to achieve the two goals: (i) 
detect enough early mergers made of building blocks not 
affected much yet by star formation and evolution and (ii) 
probe most of the GIB source population at large wave- 
lengths, we will have to do extragalactic surveys with ALMA 
using a substantial fraction of the time to find the sources. 
Such large surveys including the whole ALMA collaboration 
would be much more efficient in terms of scientific progress 
than smaller area surveys conducted by individual smaller 
teams. 



8 SUMMARY 

We have developped a phenomcnological model that con- 
strains in a simple way the IR luminosity function evolution 
with the redshift, and fits all the existing source counts and 
redshift distribution, GIB intensity and for the first time GIB 
fiuctuations observations from the mid-IR to the submm 
range. The model has been used to give some predictions 
for future Herschel deep survey observations and the all- 
sky Planck survey. It comes out that the planned experi- 
ments {SIRTF, Herschel, Planck) will be mostly limited by 
the confusion. To find out a large number of objects that 
dominate the LF at high redshift (z>2), future experiments 
need both the angular resolution and sensitivity. This can 
be achieved in the submm only thanks to interferometres 
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such as ALMA. However, mapping large fractions of the sky 
with high signal-to-noise ratio will take a lot of time (for ex- 
ample, 1.9 years to map 5 square degrees that resolve 50% 
of the CIB at 1.3 mm with ALMA). 



APPENDIX A: 

We provide, in an electronic form through a web page**, 
a distribution of the model's outputs and programs (to be 
used in IDL) containing: 

• The array dN/(dlnL.dz) as a function of L and z, dS/dz 
as a function of L and z and Si/ in Jy for each luminosity 
and redshift, from 10 to 2000 /jm and the evolution of the 
LF for both the normal and the starburst populations (for 
fiA=0.7, f7(,=0.3 and h=0.65) 

• Some uscfull programs that compute the integral 
counts, detected sources, CIB and fluctuations redshift dis- 
tribution, and the level of the fluctuations from the previous 
data cubes. 
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Table 2. Predicted CIB intensity at 15 ^lnl (ISOCAM filter), 60 and 100 ftm (IRAS filters), 170 /im {ISOPHOT filter), 350 ^tm. (filter 
such as AA / A= 1/3) and 850 fj,m {SCUBA filter) compared to measurements. 



A 

;um 


Predicted CIB 
MJy/sr 


Predicted CIB 
Wm-^sr-^ 


Measured CIB 
Wm-^sr-i 




15 


1.25 10-2 


2.5 10-8 


> 2.4±0.5 10-9 


(1) 


60 


0.12 


6.1 10-9 






100 


0.35 


1.1 10-* 


~1.5 10-8 


(2) 


170 


0.76 


1.3 10-® 


1.4±0.3 10-8 


(3) 


350 


0.76 


6.5 10-9 


5.63l^ j° 10-9 


(4) 


850 


0.20 


6.9 10-10 


5.041^-31 10-10 


(4) 



(1) Elbaz et al. (2002) 
(2) Renault et al. (2001) 
(3) From Kiss et al. (2001) and extrapolation of DIRBE measurements 
(4) Prom Fixsen et al. (1998) and Lagache et al. (1999) 



Table 3. Comparison of the predicted CIB fluctuations (for S < Smax m Jy^/sr) and the observations. 



A (/um) 


Smax (mjy) 


Observations (Jy^/sr) 


References 


Model (Jy^/sr) 


170 


1000 


~25000 


Sorel et al., in prep 


23694 


170 


250 


13000±3000 


Matsuhara et al. 2000 


15644 


170 


100 


7400 


Lagache & Puget 2000 


11629 


100 


700* 


5800±1000 


Miville-Deschcnes et al. 2002 


10307 


90 


150 


12000±2000 


Matsuhara et al. 2000 


5290 


60 


1000 


1600±300 


Miville-Deschenes et al. 2002 


2507 



* Bright sources in Miville-Deschenes et al. (2002) are removed at 60 and 100 iJ,m using the cut of 1 Jy at 60 /im. Since the 60 /im 
sources are mostly starburst galaxies, we estimate that 1 Jy at 60 /um is equivalent to 0.7 Jy at 100 /um. 



Table 4. PACS and SPIRE la confusion noise and its associated flux limit (Sum) for the photometric criterion (qp^ot=5) and the source 
density criterion (with the equivalent 'idensity)- 









(7 




S;im 








(mJy) 




(mJy) 


PACS 7b ixm. 


(Iphot 


=5.0 


2.26 10" 


3 


1.12 10-2 




^density ~ 


= 8.9 


1.42 10" 


2 


1.26 10-1 


PACS 110 ijoa 


(Iphot 


=5.0 


1.98 10" 


2 


1.0 10-1 




^density~ 


= 8.7 


1.02 10" 


-1 


8.91 10-1 


PACS 170 ixm. 


(Iphot 


=5.0 


3.97 10" 


-1 


2.00 




^density 


7.13 


9.93 10- 


-1 


7.08 


SPIRE 250 ixm 


(Iphot 


=5.0 


2.51 




12.6 




^density 


=5.2 


2.70 




14.1 


SPIRE 350 ixm 


(Iphot 


=5.0 


4.4 




22.4 




^density ~ 


= 3.6 


3.52 




12.6 


SPIRE 550 urn 


Iphot 


=5.0 


3.69 




17.8 




^density 


=2.5 


3.18 




7.94 
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Table 5. Designed surveys that could be done with SPIRE (Numbers are from the 350 /im channel). 



Surface 
(Sq. Deg.) 


(mJy) 


(mJy) 


5<Ttot 

(mJy) 


Days 


Number of 
sources 


% resolved 
CIB 


400 


100 


28.21 


103.9 


18 


4768 


1 


100 


15.3 


22.4 


27.1 


192 


33451 


6.7 


8 


7.5 


22.4 


23.6 


64 


3533 


7.8 



Unresolved sources below 5<Ti„st= 100 mJy induce a confusion noise of iTcon/=5.63 mJy. 



Table 6. Designed surveys that could be done with PACS. 



Surface 


A 

(/im) 


Days"^ 


(mJy) 


q . 6 
(mjy) 


Number of 
sources 


% resolved 
CIB 


20 Sq. Deg. 


170 


88 


7.08 


10.01 


87 322 


48.7 


625 Sq. Arcmin 


110 


67 


0.89 


1.26 


1955 


77 


25 Sq. Arcmin 


75 


96 


0.13 


0.18 


192 


87 



" Depending on the scanning/ chopping/beam switching strategy, there may be some overhead of about 20% 



Table 7. Planck sensitivities (5(Ti„st), confusion limit {Siim= 5aconf)} confusion induced by sources between S(j^ and bcTinst {^cTadd)^ 
and total 5o" noise (which is the quadratic sum of the contribution from the detector noise and the extragalactic confusion noise). Also 
given are the cold, starburst and total galaxy densities. 



A 
(/im) 


(mJy) 


(mJy) 


(mJy) 


5(Ttot 

(mJy) 


Nco(d(S> batot) 
(/sr) 


Nss(S> 5atot) 
(/sr) 


N(S> 5<7tot) 
(/sr) 


350 


216.5 


447 





497 


1342 


40 


1382 


550 


219 


200 


7.9 


297 


187 


15 


202 


850 


97 


79.4 


3.2 


125 


72 


8 


80 


i:-!8() 


■ )~ .■) 


22.1 


2.() 


02 


.">) 


1 


:-!<) 


2097 


41.5 


11.2 


2.4 


43 


23 


3 


26 
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Table 8. Sky area in square degrees to be covered to detect more than 100 sources in a redshift range Az/z=0.3, for 3 lO" Lq (top) 
and 3 lO^'^ Lq (bottom) starburst galaxies. 



3 10" Lq 


z 


S350 (mjy) 


Ssso (mJy) 


S1300 (mJy) 


Surface 




1 


2.95 


0.30 


0.09 


0.7 




3 


0.78 


0.25 


0.08 


1.8 




5 


0.21 


0.26 


0.11 


10 




7 


0.065 


0.21 


0.12 


60 




3 10^2 Lq 


z 


S3BO (mJy) 


Ssso (mJy) 


S1300 (mJy) 


Surface 




1 


22.1 


2.05 


0.7 


5.7 




3 


7.51 


1.82 


0.6 


6.1 




5 


2.73 


2.07 


0.8 


53 




7 


1.01 


1.89 


0.9 


398 



Table 9. Angular resolution and telescope diameter needed to have the confusion limit at a flux level Smin such that sources above 
Smin contribute for about 10, 30, 60 and 80 % of the CIB (top: 350 /xm, middle: 850 /xm, bottom; 1300 /im). 



% CIB 


S350 


(mJy) 


Log N (/sr) 


6° (arcsec) 


D*" (metres) 


80 


0.9 




8.01 


3.2 


23 


60 


2.5 




7.71 


4.8 


15.1 


30 


8 




7.09 


10 


7.2 


10 


18 




6.32 


23.7 


3.0 




% CIB 


S850 


(mJy) 


Log N (/sr) 


6° (arcsec) 


D*" (metres) 


80 


0.05 




8.35 


1.6 


113 


60 


0.2 




8.03 


3.0 


59 


30 


1 




7.45 


7.0 


25 


10 


2 




6.77 


15.3 


12 




% CIB 


Siaoo 


(mJy) 


Log N (/sr) 


(arcsec) 


D*" (metres) 


80 


0.02 




8.38 


1.6 


173 


60 


0.06 




8.10 


2.7 


101 


30 


0.3 




7.48 


6.6 


40 


10 


0.7 




6.75 


15.6 


17 



" 9 has been computed using NxS^ = 1/30 
D = A/ 6> 
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Figure 7. Co-moving evolution of tiie luminosity function. The dotted line is for the normal galaxies and the dotted-dashed line, for 
the starburst galaxies. The continuous lines corresponds to both starburst and normal galaxies and the dashed line is the LF at z=0 for 
compajrison. 
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Redshift 

Figure 8. Co-inoviug lurriiiiosity density distribution for the starburst galaxies (dash line), normal galaxies (dot-dash line), and both 
normal and starburst galaxies (continuous line). Also shown for comparison is the co- moving luminosity density distribution from all 
cases of Gispert et al. 2001 (crosses with error bars), together with the best fit passing through all cases (dot line). 
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Figure 9. Number counts at 15, 60, 170 and 850 ;tni (in Log) together with the model predictions (starburst galaxies: dash line, normal 
galaxies: dot-dash line, both normal and starburst galaxies: continuous line). Data at 15 /xm are from Elbaz et al. (1999), at 170 /xm from 
Dole et al. (2001), at 60 ixm from Hacking & Houck (1987), Gregorich et al. (1995), Bertin et al. (1997), Lonsdale et al. (1990), Saunders 
et al. (1990) and Rowan-Robinson et al. (1990) and at 850 /im from Small et al. (1997), Hughes et al. (1998), Barger et al. (1999), Blain 
et al. (1999), Borys et al. (2001), Scott et al. (2002) and Webb et al. (2002). 
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Figure 10. Predicted redshift distribution for resolved sources at 15, 60, 170 and 850 /im. 
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Figure 11. Cosmic Background from mid-IR to the millimetre wavelength. The 6.5 (Desert, private communication), 12 (Clements ct al. 
1999) and 15 jum (Elbaz et al 2002) lower Umits come from ISOCAM number counts; the upper limit "CAT" is from Renault et al. (2001) 
and the cross upper limits W/B are from Biller et al. (1998). At longer wavelength, we have the 60 /im estimate from Miville-Deschenes 
et al. 2002 (A), the upper limit from Finkbeiner et al. (2000), and the lower limit from number counts at 60 /im (Lonsdale et al. 1990); 
at 140 and 240 fxm are displayed the Lagache et al. 2000 (o) and Hauser et al. 1998 (*) DIRBE values; at 100 /xm is given the lower 
limit from Dwck ct al. (1998) together with the estimates from Renault ct al. 2001 (+) and the determination of Lagache ct al. 2000 (o); 
at 170 (Puget et al. 1999) and 850 /im (Barger et al. 1999) are lower limits from number counts. The analytic form of the CIB at the 
FIRAS wavelengths is from Fixsen et al. (1998). The CIB derived from the model at selected wavelengths is given by filled squares. 
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Figure 12. Redshift distribution of sources making the CIB intensity and its fluctuations (S<135 mJy) at 170 jum. 
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Figure 14. SPIRE 350 //m redshift distribution for sources detected in the very-large area survey (continuous line, S>103.9 mjy, 
multipUed by 2) and the confusion-Umited survey (dashed line, S>27.1 mJy). 
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Figure 15. Number of starburst galaxies that can be detected in different redsfiift ranges (with a Az/z=0.5) at 350 fim by the confusion 
hmited survey of 100 square degrees as a function of bolonietric luminosity. The horizontal line shows the 50 sources needed in a Az/z=0.5 
bin for ~14% accuracy. The plots are limited to the fluxes above the detection limit and to luminosities below 10^'* Lq. 
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Figure 16. Prom bottom to top: redshift distribution for the ultra-deep, deep and shallow PACS surveys at 75, 110 and 170 /um 
respectively. 
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